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ABSTRACT 

Q . Whereas a nonrelativistic distorted wave model fails to quantitatively describe analyzing 

CSJ ! power data for exclusive proton- induced proton-knockout from the 3si/2 state in ^°^Pb at 

202 MeV, the corresponding relativistic prediction provides a perfect description, thus sug- 
! gesting that the Dirac equation is the more appropriate underlying dynamical equation. We 
check the consistency of this result by comparing predictions for both dynamical models to 
new high resolution data for 3si/2 knockout in ^osp]^ ^ higher incident energy of 392 MeV. 
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(3 . Recently we demonstrated that a relativistic distorted wave model, represented by the 

solid line in Fig. (1), provides a perfect description of the energy-sharing analyzing power 
Q ■ data for exclusive proton knockout from the 3si/2 state in ^°^Pb at an incident energy 
P5 , of 202 MeV. The corresponding nonrelativistic prediction, indicated by the dashed line in 
Fig. (1), completely fails to describe the data despite exhaustive corrections including differ- 
ent kinematic prescriptions for the NN amplitudes, non-local corrections to the scattering 
wave functions, and density-dependent modifications to the free NN scattering amplitudes 
[1, 2]. On the other hand, both Dirac and Schrodinger-equation-based models provide an 
excellent description of the unpolarized energy-sharing cross section, thus highlighting the 
application of spin observables, such as the analyzing power, for discriminating between 
different dynamical effects in nuclear systems. 

Before claiming with absolute certainty that the relativistic Dirac equation is indeed the 
appropriate underlying dynamical equation, it is essential to check the consistency of the 
202 MeV result by considering 3si/2 knockout at another incident energy. In particular, 
one needs to evaluate to which extent the relatively poor energy resolution (310 KeV at 
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FIG. 1: Energy-sharing analyzing power for the knockout of protons from the 3si/2 state in 
^°^Pb, at an incident energy of 202 MeV, for coincident coplanar scattering angles {9a' = 28.0°, 
61, = —54.6°), plotted as a function of the kinetic energy of the proton scattered at angle 6a'- The 
dashed and solid curves represent the nonrelativistic and relativistic distorted wave predictions 
respectively. The data are from Ref. [1]. 

FWHM), and resulting offline peak fitting analysis associated with the 202 MeV data, in- 
fluences the quality of 3si/2 analyzing power data and interpretation thereof [1]. Indeed, 
the Kyushu University experimental nuclear physics group have recently performed a high 
resolution (250 KeV at FWHM) study of the ^^^Ph{p, 2pf^'^Te reaction at an incident en- 
ergy of 392 MeV at the Research Center for Nuclear Physics in Japan [3] : cross section and 
analyzing power data (energy-sharing and angular distributions) were measured for proton 
knockout from the 3s 1/2- , 2^3/2, 2^5/2 and l/iii/2-states. A typical binding-energy spectrum 
is displayed in Fig. (2). 

In this paper, we focus on the theoretical interpretation of the 3si/2 cross section and 
analyzing power data at 392 MeV. In particular we compare both relativistic and nonrela- 
tivistic distorted wave predictions to these new data so as to confirm or refute the previous 
claim at 202 MeV that relativistic dynamics are essential for a correct description of the 3si/2 
analyzing power. Furthermore, we also identify additional independent polarization transfer 
observables which need to be measured to further study the role of relativity in nuclear re- 
actions. Note that the original motivation for choosing a heavy target nucleus such as ^^'^Pb 
was to maximize the influence of the nuclear medium of the scattering wave functions, while 
still maintaining the validity of the impulse approximation, and also avoiding complications 
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FIG. 2: A typical binding-energy spectrum for the 2°^Pb(p, 2pf°'^Te reaction at 392 MeV. The 
relevant single-particle proton-hole states in ^'^^Pb are indicated. The data are from Ref. [3]. 

associated with the inclusion of recoil corrections in the relativistic Dirac equation. In ad- 
dition, for the knockout of 3si/2 valence protons one expects density-dependent corrections 
to the NN interaction to be negligible. 

This paper is organized as follows: we briefly review the main aspects of the relativistic 
and nonrelativistic distorted wave models in Sees. (II) and (III) respectively. The relevant 
scattering observables are defined in Sec. (IV), and results and conclusions are presented in 
Sec. (V). 

II. RELATIVISTIC MODEL 

The formalism for the relativistic distorted wave model has been presented in Refs. [2, 4- 
7] . In this paper we briefly allude to the most important aspects of the model and refer the 
interested reader to the latter references for more detail. 

The exclusive (p, 2p) reaction of interest is schematically depicted in Fig. (3), whereby an 
incident proton, a, knocks out a bound proton, b, from a specific orbital in the target nucleus 
A. resulting in three particles in the final state, namely the recoil residual nucleus, C, and 
two outgoing protons, a' and b, which are detected in coincidence at coplanar laboratory 
scattering angles, 9a' and 9b, respectively. All kinematic quantities are completely determined 
by specifying the rest masses, m^, of particles, where i — {a,A, a', b, C), the laboratory kinetic 
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FIG. 3: Schematic representation for the coplanar (p, 2p) reaction of interest. 

energy, T^, of incident particle a, the laboratory kinetic energy, T^/, of scattered particle a', 
the laboratory scattering angles 6a' and 6^, and also the excitation energy of the proton that 
is to be knocked out of the target nucleus, A. 

We adopt a zero-range approximation to the NN interaction, whereby the relativistic 
distorted wave transition matrix element is given by 

TLJMj{Sa,Sa',Sb) = j df[^'^~\f,kaiC,Sa')®^^'\r,hc,Sb)] 

FM^i^, ^cT) [ Ka, Sa) ® $f,M,(r ) ] (1) 

where denotes a Kronecker product. The four-component scattering wave functions, 
ip{r, ki, Si), are solutions to the fixed-energy Dirac equation with spherical scalar and time- 
like vector nuclear optical potentials: '^^~^\r, kaA, Sa) is the relativistic scattering wave func- 
tion of the incident particle, a, with outgoing boundary conditions [indicated by the su- 
perscript (+)], where kaA is the momentum of particle a in the (a + A) center-of-mass 
system, and Sa is the spin projection of particle a with respect to kaA as the ^-quantization 
axis; ^^ \f,kjc,Sj) is the adjoint relativistic scattering wave function for particle j [ j = 
{a',b)] with incoming boundary conditions [indicated by the superscript (— )], where kjc 
is the momentum of particle j in the (j + C) center-of-mass system, and sj is the spin 
projection of particle j with respect to kjc as the z-quantization axis. We employ a global 
scalar and vector Dirac optical potential parameter set for the distorting optical potentials. 
More specifically, we use the energy-dependent mass-independent "EDAI-fit" parameter set 
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which has been constrained by proton elastic scattering data on ^'^^Pb for incident proton 
energies between 21 MeV and 1040 MeV [8]. The four-component relativistic boundstate 
proton wave function, which is denoted by $f jmj(^) Eq. (1) and labeled by single-particle 
quantum numbers L, J, and Mj, is obtained via selfconsistent solution to the Dirac-Hartree 
field equations within the context of the relativistic mean field approximation associated 
with the QHDII Lagrangian density [9] of quantum hadrodynamics [10]. We also employ 
the impulse approximation which assumes that the form of the NN scattering matrix in 
the nuclear medium is the same as that for free NN scattering. In addition, we adopt the 
lAl representation [11] which parameterizes the NN scattering matrix Fnn{T^^^,61'^) in 
terms of five Lorcntz invariant amplitudes (scalar, pseudoscalar, vector, axial-vector, ten- 
sor) which are directly related to the five nonrclativistic Wolfenstein amplitudes describing 
on-shell NN experimental scattering data: T^q° and O^q represent the effective two-body 
laboratory kinetic energy and center-of-mass scattering angles, respectively, based on the 
so-called final-energy prescription [12]. 

III. NONRELATIVISTIC MODEL 

All of the nonrclativistic calculations are based on the computer code THREEDEE of 
Chant and Roos [12], where the nonrclativistic transition amplitude, based on a zero-range 
approximation to the NN interaction, is 



where the ■^'s represent the appropriate incoming and outgoing nonrclativistic two- 
component scattering wave functions and (p is the nonrclativistic wave function of the bound 
proton to be knocked out, and ^ — A/{A+1), where A being the target mass number: both 
scattering and boundstate wave functions are solutions to the Schrodinger equation. For con- 
sistency between the relativistic and non-rclativistic calculations, the nonrclativistic radial 
boundstate wave function is approximated by taking the upper component radial wave func- 
tion of the relativistic four component boundstate wave function employed in the relativistic 
predictions. The scattering wave functions are solutions to the nonrclativistic Schrodinger 
equation employinga Schrodinger-equivalent representation [8] of the relativistic global opti- 
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cal potentials mentioned in Sec. (II). I n this paper we employ experimental amplitudes for 
the NN scattering matrix tNN which are determined from the NN Arndt phase shift anal- 
ysis (January 1999) [13]: these amplitudes are directly related to the relativistic Lorentz 
invariant amplitudes [11]. 



IV. SCATTERING OBSERVABLES 



The spin observables of interest are denoted by Dj/j and are related to the probability 
that an incident beam of particles a with spin-polarization j induces a spin-polarization i' 
for the scattered beam of particles a': the subscript j — (0, £, n, s) is used to specify the 
polarization of the incident beam, a, along any of the orthogonal directions 

i ^ z ^ kaA, n ^ y ^ kaAX ka'c , s ^ X ^ n X £ , (3) 

and the subscript i' = (0, n', s') denotes the polarization of the scattered beam, a', along 
any of the orthogonal directions: 

i' = = k^,Q ^ fi' ^ h ^ y , s' ^ x' ^ h X £' . (4) 

The choice j {i') = is used to denote an unpolarized incident (scattered) beam. With the 
above coordinate axes in the initial and final channels, the spin observables, Dj^j, are defined 

by 

^ _ Em,,., Tr(Ta,Tta,0 

where D„o = P refers to the induced polarization, Do„ = A^^ denotes the analyzing power, 
and the other polarization transfer observables of interest are D„„, D^/^, D^/^, D^/^, and D^/^. 
The denominator of Eq. (5) is related to the unpolarized triple differential cross section, i.e., 

^ ^ JO ^ SijCTcalc, (6) 

ala' aiZa' OlZft 

where Fkin is a kinematic factor and Slj is the spectroscopic factor [5, 6]. In Eq. (5), the 
symbols Uj/ and aj denote the usual 2x2 Pauli spin matrices, and the 2x2 matrix T is 
given by 

T = I ^'-+1 . ,-1 ^-i . I (8) 



where = i| and s^/ = 



±1 refer to the spin projections of particles a and a' along the 



z and z' axes, defined in Eqs. (3) and (4), respectively; the matrix T^y^"' is related to the 
relativistic {p,2p) transition matrix element TLjMj{sa, Sa', Sb), defined in Eq. (1) via 



V. RESULTS AND CONCLUSIONS 

The experimental data (experiment E205) were measured using the dual-arm spectrome- 
ter at RCNP [3]. For a direct comparison to data both relativistic (RDWIA)and nonrelativis- 
tic (NRDWIA) distorted wave predictions were corrected (using a Monte-Carlo simulation) 
for the finite angle acceptance of both Grand Raiden and the Large Acceptance Spectrom- 
eters: note that the current version of the relativistic code does not consider out-of-plane 
predictions, and consequently the RDWIA calculations exclude finite azimuthal angle 
corrections to the sohd angle. 

We now compare theoretical RDWIA and NRDWIA energy-sharing and angular distri- 
bution predictions of cross sections a and analyzing powers Aj^ to experimental data for 
3si/2 knockout from ^osp^^ MeV. In Fig. (4) the energy-sharing distributions (left 

panel) are plotted as a function of the kinetic energy Ta' for coincident coplanar laboratory 
scattering angles (32.5°, —50.0°), and the angular distributions (right panel) are plotted as 
a function of the scattering angle 9b, for T^/ = 250 MeV and 6'„/ = 32.5°. Both RDWIA 
(solid curves) and NRDWIA (dashed curves) models provide a satisfactory description of the 
shape of the unpolarized cross sections for both energy- sharing and angular distributions. 
The analysis regarding the extraction of spectroscopic factors Slj and corresponding error 
bars is currently in progress: the values of Slj represent single-particle state occupation 
numbers and are obtained by normalizing the calculated cross sections (Jcaic [see Eq. (7)] to 
the experimental cross section data. 

Next we turn our attention to the analyzing power. As is the case for the 202 MeV 
data, the RDWIA energy-sharing distribution, indicated by the solid line in Fig. (4), is 
significantly reduced compared to the corresponding NRDWIA calculations (dashed line). 
Furthermore, the RDWIA prediction provides a better overall quantitative description of 
the data. Recall that for the 202 MeV data the energy-sharing analyzing power distribution 
for RDWIA was consistently reduced compared to the RDWIA calculations for all values of 




TLJMj{Sa, Sa', Sb) ■ 
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Ta'. The latter is also true at 392 MeV, except at the maximum value of the cross section 
(at Ta' ^ 245 MeV) where both RDWIA and NRDWIA describe the data equally well. The 
RDWIA angular distribution is also consistently reduced compared to the NRDWIA result, 
except at ~ 50° (corresponding to the maximum value of the cross section) where both 
models descrribe the data equally well. Hence, in general we conclude that at both 202 and 
392 MeV, the RDWIA model is superior compared to the NRDWIA model, thus suggesting 
that relativistic dynamics are important for for describing (p, 2p) reactions. 




T,. (MeV) 6^ (deg) 

FIG. 4: Unpolarized triple differential cross section a and analyzing power Ay for proton knockout 
from the 3si/2 state in ^'^^Pb at an incident energy of 392 MeV plotted as a function of the kinetic 
energy T^/ for coincident coplanar laboratory scattering angles (32.5°, —50.0°) (left panel), and 
as a function of the scattering angle Ob, for Ta' = 250 MeV and 6a' = 32.5° (right panel). The 
dashed and solid curves represent the nonrelativistic and relativistic distorted wave predictions 
respectively. The data are from Ref. [3]. 

It is desirable to check the consistency of this result by also measuring other polariza- 
tion observables which are sensitive to differences between relativistic versus nonrelativistic 
dynamical models. Based on the same kinematic conditions as for the 392 MeV analyzing 
power, we also identify the spin observables, D„„, Dg'g and D^/^ as good candidates for fur- 
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ther studying the role of relativity in nuclei: see Fig. (5): Other spin observables are not 
displayed since they are less sensitive to different dynamical effect. 
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FIG. 5: Polarization transfer observables, D„„, Dg's and De'e, plotted as a function of either T^' 
(left panel) or Oi, (right panel) for kinematics corresponding to observables plotted in the left and 
right panels of Fig. (4) respectively. The dashed and solid curves represent the nonrelativistic and 
relativistic distorted wave predictions respectively. 



The success of the RDWIA model to describe the knockout of the 3si/2 valence nucleons in 
208p|-j inspires confidence to extend our model to systematically address the topical question 
of how the free NN interaction is modified by the presence of the neighbouring nucleons in 
nuclei: the exclusive nature of {p, 2p) reactions allows one to selectively knockout protons 
from deep- to low-lying singlc-particlc orbitals in nuclei, thus enabling one to systematically 
study the effect of the nuclear density on the NN interaction. In addition, we intend to 
apply our RDWIA model to study exclusive proton knockout reactions from proton-rich 
and neutron-rich exotic nuclei using inverse kinematics at future radioactive beam facilities 
such as RIKEN and GSI. A drawback of the current implementation of the RDWIA model 
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is the use of the ambiguous lAl parameterization for the NN scattering matrix. Future work 
will include study the effect of replacing the lAl with the unambiguous IA2 representation in 
terms of 44 independent invariant amplitudes (of which the lAl representation is a subset) 
which are consistent with parity and time-reversal invariance as well as charge 
symmetry [14]. 
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